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Abstract 

We study for the first time the complete sub-millimeter spectra (450 GHz to 1550 GHz) of a 
sample of nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer 
(SPIRE/FTS) onboard Herschel". The GO ladder (from Jup = 4 to 12) is the most promi- 
nent spectral feature in this range. These CO lines probe warm molecular gas that can be 
heated by ultraviolet photons, shocks, or X-rays originated in the active galactic nucleus 
or in young star-forming regions. In these proceedings we investigate the physical origin 
of the CO emission using the averaged CO spectral line energy distribution (SLED) of six 
Seyfert galaxies. We use a radiative transfer model assuming an isothermal homogeneous 
medium to estimate the molecular gas conditions. We also compare this CO SLED with 
the predictions of photon and X-ray dominated region (PDR and XDR) models. 



"'Herschel is an ESA space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from NASA. 
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1 Introduction 

To date, the sub-millimeter (sub-mm) spectrum of galaxies has not been fully explored, 
in part due to the low transparency of the atmosphere in this range. The SPIRE Fourier 
transform spectrometer (SPIRE/FTS) [2, 6, 12] onboard the Herschel Space Observatory [8] 
fills this gap, covering the spectral range between 450 and 1440 GHz (210 and 670 fim). 

In the SPIRE/FTS spectra the most obvious spectral feature is the mid-J CO ladder 
(<^up = 4 to 13). CO is one of the most abundant molecules in the interstellar medium; 
therefore it is considered a good tracer of the molecular gas phase. Specifically, the mid-J 
CO line emission originate from warm molecular gas (upper-level energies ranging from 55 
to 500 K above the ground state) that can be excited by ultraviolet (UV) photons in photon 
dominated regions (PDRs; e.g., [14]), X-rays in X-ray dominated regions (XDRs; e.g., [5]) or 
by shocks (e.g., [1]). 

In these proceedings we analyze the average CO spectral line energy distribution 
(SLED) of local Seyfert galaxies, including both low- and mid-J CO transitions, in order 
to understand the physical conditions of the warm molecular gas and the dominant excita- 
tion mechanism in this class of galaxies. 

2 Sub-mm Spectra of Local Seyfert Galaxies 

We obtained sub-mm SPIRE/FTS spectra of eleven nearby {d = 12 — 160 Mpc) Seyfert 
galaxies selected from the 12 /im Galaxy Sample [10] through two guaranteed time projects: 
"Bright Seyfert Nuclei: FTS spectroscopy" (PI: L. Spinoglio) and "Physical Processes in the 
Interstellar Medium of Very Nearby Galaxies" (PI: C. D. Wilson). Our sample includes 
galaxies with different nuclear activity classes in similar proportions. Their infrared (IR) 
luminosities^ range from 10^^ to 10^^ Lq and their hard X-ray luminosities^ from 10"^" to 
10^3 erg s-^ The SPIRE/FTS beam size is 20-40" that corresponds to 1.4-32 kpc at the 
distance of our galaxies. More details on the sample and data reduction are given in [7]. 

In the SPIRE/FTS spectra of these galaxies we detected transitions of several molecules. 
The brightest lines are those produced by CO, but lines from other molecules like H2O, 
OH"'", and HE are often detected. In addition, atomic lines from neutral gas ([C I] at 492 
and 809 GHz) and ionized gas ([N II] at 1461GHz) are bright in these galaxies [7]. In these 
proceedings we focus on the mid-J CO emission. 

For six of these galaxies we detected more than three CO lines at >3a level in their 
SPIRE/FTS spectra (see Figure 1) that were used, together with ground observations of 
the low- J CO lines {Jnp = 1 to 3), to construct the CO SLEDs for each galaxy. The CO 
SLEDs of five of them (UGC 05101, NGC 3227, NGC 4388, NGC 7130, and NGC 7582) 
are remarkably similar. Therefore we averaged them (normalized by their CO luminosity) 
to construct an average SLED which we used to study the typical conditions in the warm 
molecular gas in Seyfert galaxies. We note that the Seyfert type, the AGN luminosity, or the 

^The IR luminosity refers to the 8-1000 fim galaxy integrated luminosity while the hard X-ray luminosity 
is the 2-10 keV luminosity. 
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Figure 1: Observed SPIRE/FTS spectra of the Seyfert 1.5 NGC 3227. The dashed red hnes 
mark the position of the detected hnes. The error bars indicate the median la uncertainty. 



total IR luminosity (mainly related to their star-formation rate) do not seem to affect the 
CO SLED shape in our sample of galaxies. 



3 CO Emission Modeling 
3.1 Radiative Transfer Model 

We used the non-LTE radiative transfer code RADEX, with the assumption that the mid- 
J CO emission originates from a homogeneous medium, to infer the physical conditions of 
the emitting regions. We produced a grid of models that was used to model the average 
CO SLED. The best fit RADEX model is shown in Figure 2 (dotted red line). The derived 
kinetic temperature is 500^250 ^ -^2 density iO'^-''='=^ '^ cm~^. As expected, these values 

agree with those obtained from the individual galaxies (see [7]). Moreover, these conditions 
are similar to those observed in the starburst M 82 [3]. Differently, in the ultra-luminous 
IR galaxy Arp 220 the gas temperature is higher (~ 1300 K) [9], whereas in the Seyfert 2 
NGC 1068 it is lower (~ 90 K) [11]. This is consistent with the mid- J CO emission from 
our sub-sample of Seyfert galaxies being driven by star-formation activity rather than by the 
AGN. 



3.2 PDR Models 

In star-forming regions where UV photons from young stars heat the ISM we expect PDR 
excitation of the molecular gas. Therefore we compared the average CO SLED with the PDR 
models of [14]. In Figure 2 we plot the best-fit model (solid green line). The plotted model cor- 
responds to = 10^'^^ cm~^ and far-UV flux 10^'^ Gq, where Go = 1-6 x 10"'^ ergcm"^ s~^. 
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Figure 2: Average CO SLED of our sample of nearby Seyfert galaxies (gray squares) nor- 
malized by the total CO luminosity. The dotted red line represents the best-fit RADEX 
model. The solid green line is the best-fit PDR model [14]. XDR models of [5] are plotted 
with purple, blue, and orange dashed lines (see Section 3.3). 



However, models with = 10^'^ — 10^'^ cm~^ and -Ffuv = 10 — 10^'^ Gq^ also provide 
equally good fits. These PDR models provide a good fit to the emission of CO lines with 
Jup < 8, however emission from higher- J transitions is clearly underpredicted. 

Figure 3 shows the Leo / -^ir'^ ratio as a function of the incident far-UV fiux and gas 
density. In our galaxies this ratio is ~10~^, however this is a lower limit because not all 
the IR luminosity have to be associated with the warm molecular gas emitting the mid-J 
CO lines. Consequently, from this figure we can constrain the far-UV radiation field to be 
between 10 and 10^ Go- 



3.3 XDR Models 

We also fitted our averaged CO SLED using the XDR models of [5]. These authors model 
three possible scenarios for XDR emission: high-density (nHa ~ 10^ — lO^'^cm"^) compact 
molecular clouds close to the AGN (XDR A); "traditional" mid-range density {nn^ ~ 10^ — 
10^ cm~'^) molecular clouds (XDR B); and diffuse molecular gas (nHa ~ 10^ — 10^ cm~'^) where 
XDR B type clouds are usually embedded (XDR C). The best-fit models are plotted in Figure 
2 with dashed lines. The XDR models better reproduce the CO emission for Jup > 8 than 
do the PDR models. The best one is the XDR B, although it systematically underpredicts 
the CO luminosities for Jup > 8. Conversely, the XDR A model underpredict CO lines with 
•Aip < 6, but agrees with the observed higher- J emission. 



■^fpuv = lO'''^ Go is the upper limit of the grid of PDR models, so higher -Fpuv could be compatible as 
well. 

^Lco is defined as the integrated luminosity of the CO lines from Jup = 4 to Jup = 12. 
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Figure 3: The black contours represent the logarithm of the Lqo /^ir ratio as a function 
of the incident far-UV flux and gas density. The models capable of fitting the observed CO 
SLED (up to Jup = 7) are marked with a shaded green area. We assumed Lir, = 2 x Lfuv [4]- 

The Lqo/Lx^ ratio, for these two XDR models is ~10~^; therefore for our Seyfert 
galaxies the observed Lqo implies X-ray luminosities 100 — 10 000 times higher than those 
observed. Other XDR models in these grids have higher Lqq/Lx ratios (up to ~10~^) that 
would be compatible with the AGN luminosities although they produce worse fits to the CO 
SLED. On the other hand we can reject the extended diffuse XDR emission (C model) since 
the observed CO surface brightness in our galaxies (>10~^ ergcm"^ s~^ sr~^) is almost two 
order of magnitude higher than that of the best-fit model (5 x 10"^ ergcm^^ s^^ sr~^). 

4 Summary and Conclusions 

We obtained the CO SLED (from Jup = 1 to 12) of six nearby Seyfert galaxies using SPIRE/ 
FTS spectroscopic observations. The average CO SLED was fitted using simple radiative 
transfer models (RADEX) as well as more elaborate PDR and XDR models. 

From the RADEX modeling we derived a kinetic temperature of 500!!^250 ^ ^^"^ -^2 
density of lO^''^^^''^ cm~'^. However, we note that an homogeneous model with a single tem- 
perature and gas density may not be a realistic representation of the molecular gas conditions 
in a Seyfert galaxy. 

PDR models can explain the observed CO SLED up to Jup ~ 8. Higher J CO lines 
are underpredicted by these models and another component (PDR or XDR) is needed to 
reproduce the SLED. We also show that XDR models can fit the CO SLED, however the 
measured Lqo/Lx ratios disagree with those expected for these models. Nevertheless, there 



'^Lx is the 1-100 keV X-ray luminosity. 
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are XDR models compatible with the observed Lqq I Lx ratios, although they do not fit well 
the CO SLED. In summary, an individual PDR or XDR model cannot explain the observed 
CO emission in these local Seyfert galaxies and, at least, two components (two PDRs or 
PDR+XDR) are needed. 

Future high sensitivity and high angular resolution observations with ALMA of the 
mid-J CO lines will help to disentangle the AGN and star-formation CO emission, as well 
as, to determine the excitation mechanism and heating source of the warm molecular gas. 
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